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ABSTRACT 



Context. The low mass protostar IRAS 16293-2422 is a prototype Class source with respect to the studies of the chemical structure 
during the initial phases of life of Solar type stars. 

Aims. In order to derive an accurate chemical structure, a precise determination of the source physical structure is required. The scope 
of the present work is the derivation of the structure of IRAS 16293-2422. 

Methods. We have re-analyzed all available continuum data (single dish and interferometric, from millimeter to MIR) to derive 
accurate density and dust temperature profiles. Using ISO observations of water, we have also reconstructed the gas temperature 
profile. 

Results. Our analysis shows that the envelope surrounding IRAS 16293-2422 is well described by the Shu "inside-out" collapsing 
envelope model or a single power-law density profile with index equal to 1.8. In contrast to some previous studies, our analysis does 
not show evidence of a large (> 800 AU in diameter) cavity. 

Conclusions. Although IRAS 16293-2422 is a multiple system composed by two or three objects, our reconstruction will be useful to 
derive the chemical structure of the large cold envelope surrounding these objects and the warm component, treated here as a single 
source, from single-dish observations of molecular emission. 
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1. Introduction 



stars: individual: IRAS 16293-2422 - ISM: molecules - ISM: abundances - stars: circumstellar matter 



Understanding how our Sun and Solar System formed is ar- 
guably one of the major goals of modern astrophysics. Many 
different approaches contribute to our understanding of the past 
history of the Solar System. Analyzing the relics of the ancient 
eons, comets and meteorites, is one. Studying present day ob- 
jects similar to what the Sun progenitor is another. Here we pur- 
suit the latter approach and analyze in detail the case of one of 
the best studied solar type protostars, IRAS 16293-2422 (here- 
inafter IRAS 16293). IRAS 16293 is a Class protostar in the p 
Ophiuchus complex at 120 pc from the Sun dLoinard et alJ2008h 
and has played the role of a prototypical solar-type protostar 
for astrochemical studies, just as Orion KL has done for high- 
mass protostars. This is beause of its proxi mity and the result- 
ing line strength of molecular emission (e.g . ^Walke r et al.|ll986 ; 
iMundv et alj|l992l: iBlake et al.|[T994t Ivari Dishoeck etal.lll995 , 



to mention just a few representative works from the previous 
decades). It is in this source that the phenomenon of the "super- 
deuteration'Qhas been first discovered, with the detection of sur- 
prising abundant multiply deuterated molecules: formaldehyde. 



hydrogen sulfide, methanol and wat er ("Ceccarelli et al.l [l998l : 
IVastel et alj I200I iParise et al.1 120031: Butneret al. 2007) It is 
also in this source that the first hot corino ( Bottinelli et al]l2004l : 
Ceccarelli et al. 2007) has been discovered, with the detection 
of several abundant complex organic molecules in the region 
whe re the du st grain mantles sublimatej Ceccarelli et al.ll2000bt 
Caz aux et alJ 2003: Bottinelli et al. 20Q3. 

Not surprisingly, therefore, IRAS 16293 has been the target 
of several studies to reconstruct its physical structure, na mely the 
dust and gas density and te mperature profiles ( Ceccare lli et al.l 
120003"; 'S choieretalJ I200I l2004t IJ0rgensen et all I200"5h . the 
mandatory first step to correctly eva luate the abundance of 
molecular species across the envelope. ICeccarelli et alJ (l2000ah 
used water and oxygen lines observations obtained with the 
Infrared Space Observatory (ISO) to derive the gas and dust den- 
sity and temperature profile. Conversely, I Schoier et al.l (l2002h 
used the dust continuum observations to derive the structure 
of the envelope. Furthermore, while Ceccarelli et al. assumed 



Send offprint requests to: N. Grimier 

' The super-deuteration refers to the exceptionally high abun- 
dance ratio of D-bearing molecules with respect to their H-bearing 
isotopologues found in low mass protostars, with observed D- 
molecule/H-molecule ratios reaching the unity (see e.g. the review in 
ICeccarelli etani2007r) . 



the semi-analytical solution by Shu & co-workers (Sh ulll977t 
Adams & Shu 1986) to fit the observational data, Schoi er et al.l 
(2002) assumed a single power law for the density distribution, 
and a posteriori verified that the Shu's solution also reproduced 
the observational data. The two methods lead to similar general 
conclusions: a) the envelope of IRAS 16293 is centrally peaked 
and with a density distribution in overall agreement with the 
inside-out collapse picture ( Shu 1977); b) there is a region, about 
300 AU in diameter, where the dust mantles sublimate (giving 
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rise to the phenomenon of the hot corino, mentioned above). 
However, the two different methods, not surprisingly, also lead 
to some notable differences. For example, the gas density dif- 
fers by about a factor 3 in the region where the ice sublimation 
is predicted to occur, which leads to differences in the derived 
abundances of several molecular species. 

Subsequent studies have built on the early ones to improve 
the deriv ation of the physical structure of IRAS 16293. First, the 
study by ISchoier et alj (|2004|) . based on interferometric OVRO 
observations, concluded that the envelope has a large cent ral cav- 
ity, about 800 AU in diameter. Then, the work by Jcirgense n et alJ 
(120051) . using new SPITZER data, concluded that the envelope 
has an even larger central cavity, about 1200 AU in diameter. 
Such a large central cavity has a great impact on the whole in- 
terpretation of the hot corino of IRAS 16293, as it predicts the 
absence of the mantle sublimation region. If the predicted cavity 
is real, the observed complex organic molecules must have an- 
other origin than grain mantle sublimation from thermally heated 
dust. In addition to rising an important point in itself, the Schoier 
and J0rgensen et al. works illustrate the paramount importance 
of correctly understanding the physical structure of the source to 
assess the chemical structure and all that follows. 

For this reason, in the present work, we have re-analyzed 
the available data on IRAS 16293 from scratch, considering, 
in additio n, the most re c ent ev aluation of the distance to this 
source bv iLoinard et al.l (l2008l) (120 pc instead of 160 pc, as 
assumed in Ceccarelli et al., Schoier et al. and J0rgensen et 
al. works). This new analysis is necessary and timely because 
of two important observational projects having IRAS 16293 as 
a target: a) the unbiased spectral survey in the 3, 2, 1 and 
0.8 mm bands just obtained at the IRAM and JCMT tele- 
scopes ("The IRAS 16293-2422 Milli meter And Sub-millimeter 
Spectral Survey'Q; ICaux et al l l20T0[ in prep.), and b) the unbi- 
ased spectral survey between 500 and 2000 GHz which will be 
obtained shortly with the heterodyne instrument HIFI aboard the 
Herschel Space Observatory (HSO) launched in May 2009 (the 
Herschel Guaranteed Time Key Program CHESS — Chemical 
Herschel Surveys of Star Forming RegionQ). The two projects, 
involving large international teams, will provide an accurate cen- 
sus of the molecular inventory of IRAS 16293, the largest ever 
obtained in a solar type protostar. To convert the observations 
into an accurate chemical composition across the IRAS 16293 
envelope, the dust and gas density and temperature profiles have 
to be determined accurately first. Deriving these profiles is the 
goal of the present article. 

We conclude this section by addressing the problem of the 
binarity of IRAS 16293 and how it fits with the analysis we 
present here. As soon as interferometric observations became 
available it w as realized that IRAS 16293 is in deed a proto- 
binary system ( IWoottenllI989l:lMundv et al.llI992h . composed by 
two sources: A (the south source) and B (the north source) sep- 
arated by ~4", i.e. -500 AU at 120 pc. While source B is the 
brightest in the continuum, source A is oft en, but not always, th e 
brightest in the molecular emission (e.g. Chandler et ai]|2005l) . 
The most recent observations show that IRAS 1 6293 is indeed a 
triple system, with the source A composed by two objects, Al 
and A2, of 0.5 and 1.5 Mq respectively (Loinard et al. 2009). 
While it is clear that the multiple nature of IRAS 16293 cannot 
be neglected in general, the two projects mentioned above in- 
volve observations with single-dish telescopes, so that much of 
the structure on small scales is smeared out in these observations. 



thttp://www-laog.obs.iijf-grenoble.fr/heberges/timasss/^ 
\http://www-taog.obs.ujf-grenoble.fr/heberges/chess/, 



Specifically, the molecular line emission will be dominated by 
the cold envelope, which fills up the telescope beam, and by any 
warm component at the interior of the envelope. The major goal 
here is to give a reliable estimate of the envelope temperature 
and density profiles of both the gas and dust components, up to 
the scales where the approximation of a spherical symmetry is 
valid. What are these scales will be discussed later on, based on 
the available observations. 

The article is organized as follows. Section |2] discusses the 
derivation of the dust density and temperature distribution, based 
on the analysis of all available continuum data. Section [3] de- 
scribes the derivation of the gas temperature profile, with the 
help of ISO data to constrain the abundance of a major gas 
coolant. Finally, Section|4]discusses and summarizes the results 
of the presented study. 



2. Dust temperature and density profiles 

2.1. The data set 

The present analysis is based on the continuum emission from 
the envelope that forms/surrounds the protostar IRAS 16293. 
Three types of observations are considered: maps of the emis- 
sion, spectral energy distribution (SED) and interferometric 
observations at 1 and 3 mm. All data used have been retrieved 
from archives, except of the map at 350 fim and the 1 and 3 mm 
interferometric data that we obtained in dedicated observations. 
Below we briefly describe the data used. 

i) Continuum emission profiles 

We used the maps of the dust continuum emission at 350 jum 
(obtained at the Caltech Submillimeter Observatory; CSO), 
and 450 and 850 fim (obtained at the James Clerk Maxwell 
Telescope; JCMT). 

The 450 and 850 /im maps have been retrieved from the 
JCMT archive (website). The beam sizes are 7.5" and 14.8" 
at 450 and 850 jim respectively. Based on the many previous 
JCMT published observations, the calibration uncertainty and 
noise levels are < 10% and 0.04 Jy beam"' at 850 ;um and < 

30% and 0.3 Jy beam"' at 450 //m, respectively. 

Observations of the 350 //m continuum emission toward 
IRAS 16293 reported here were carried out in 2003 February 
using the SHARC II facility bolometer camera of the Caltech 
Submillimeter Obs ervatory (CSO) on Mauna Kea in Hawaii 
(iDoweU et al.ll2003 ). SHARC II is a 12 x 32 pixel fifled array 
with a field of view of 2.6 x I.O'. The data were taken during 
excellent submillimeter weather conditions (a 225 GHz zenith 
opacity of 0.04, corresponding to less than 1 mm of precipitable 
water). The observations were carried out using the "box-scan" 
scanning mode (see http://www.submm.caltech.edu/sharc/). Five 
10 min sc ans were redu ced together using the CRUSH software 
package (iKovacsl [2008h to produce the final calibrated image. 
Telescope pointing was checked on Juno, Vesta, or Mars, imme- 
diately before or after the science observations and the measured 
offsets were applied during data reduction. The data were taken 
before the CSO Dish Surface Optimization system (DSOS) be- 
come operational. The shape of the telescope beam was deter- 
mined from the pointing images of Vesta and Mars. It contains 
a diffraction limited main beam with a FWHM diameter of 9" 
and an error beam with a FWHM diameter of 22", with relative 
peak intensities of 0.8 and 0.2, respectively. This size of the error 
beam is consistent with earlier 3 50 //m measurements using the 
SHARC I camera (lHuntedll997h . 
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Fig. 1. Continuum emission profiles at 850 jum (left bottom panel), 450 fim (right upper panel) and 350 /im (left upper panel), plus 
the SED (right bottom panel). The curves show the emission predicted by two models with different density profiles (see text for 
more details): the solid thick line is the best fit of the maps and SED data obtained by a density profile with a 2-indexes power law 
with Q'o«r=2 and q',„=1.5 in the outer and inner envelope respectively (Shu-like density distribution); the dashed line is the best fit 
obtained by a single power-law index, q'=1.8. The dotted-dashed lines show the telescope beam profiles adopted. 



We obtained continuum emission profiles as function of 
the distance from the center of the envelope, by averaging the 
map emission over annuli at the same distance. In each case 
we adopted a radial sampling corresponding to the half of the 
resolution of the instrument. The uncertainties of continuum 
emission profiles were evaluated taking into account the cali- 
bration uncertainty, noise levels and the non sphericity of the 
source. The resulting profiles, normalized to the peak emission, 
are shown in Fig. [1] Note that the 450 and 850 /um profiles are 
identical to those reported by Sc hoier et al.. (.2002) . 

ii) Spectral Energy Distribution 

Table [T] reports the SED obtained considering all the data 
available in the literature (plus the 350 yum point obtained by us; 
see above). The millimeter and sub-millimeter data points have 
been obtained by integrating the maps over the entire envelope. 
The in tegrated flux at 1.3 mm was that quoted in lSaraceno et alj 
(119961) . We reti-ieved the IRAS fluxes from flie IRAS 
Point Source Catalog v2.1 {http://irsa.ipac.caltech.edu/cgi- 
bin/Gator/nph-dd?catalog—iraspsc). 

IRAS 16293 was observed with the InfraRed Spectrograph 
(IRS) installed aboard the Spitzer Space Telescope as part 
of the "Fro m Molecular Cores to Planet Forming Disks" 
(lEvans etalJl2003t lEvans & c2d TeamI 120051) Le gacy Program 
(AOR: 11826944, PI: Neal Evans). We use the observations 
obtained with the Long-High (LH) module (20-37 pm, R - 
600) on 2004 July 29, in staring mode. The data reduction was 
performed using the c2d pipeline S15.3.0 (Lahuisetal. 2006) 
with the pre-reduced (BCD) data. Since the MIPS map at 24 pm 
shows that the IRAS 16293 emitting region (~ 30-40") is larger 
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Table 1. The Spectral Energy Distribution of IRAS 16293-2422. 

Integrated flux in Jy. Uncertainties in Jy considering the 
calibration uncertainty, the noise levels and the uncertainty on 
the source size. * Main beam of the instrument in arcsec. 

References. (1) This paper; (2) Spitzer catalog; (3) IRAS Point Source 
Catalog v2.1; (4) ,Saraceno et am996) . 



than the LH module field of view (11.1" x 22.3"), we adopted 
the full aperture extraction method in the pipeline. In addition, 
we corrected the derived flux level for the missing flux by 
comparing the IRS spectrum integrated on the MIPS bandwith 
and the integrated MIPS flux. Note that this method assumes 
a similar distribution of the source emission in the whole IRS 
wavelength interval, ~ 21-37 pm. The correction factor derived 
by this method amounts to 22 % in the flux. 

Hi) Interferometric continuum data 
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Observations of the 3 and 1.3 mm continuum were obtained 
with the IRAM Plateau de Bure Interferometer (PdBI) in 2004 
and are described in detail in Bottinelli et al. (2004). They were 
obtained in B and C configurations of the PdBI array, resulting 
in a spatial resolution of about 0.8". 

2.2. Adopted model 

We obtained the best fit to the continuum data, described in the 
prev ious section, by using the ID radiative transfer code DUSTY 
dlvez ic & Elitzur 1997), which has been extensively used in sim- 
ilar \yorks (in cluding J0rge nsen et al.ll2002l: ISchoier et al.ll2002i 
l2OO4tlj0rgensen et al..,2005i) . Briefly, giving as input the temper- 
ature of the central object and a dust density profile, DUSTY 
self-consistently computes the dust temperature profile and the 
dust emission. A comparison between the computed 350, 450, 
850 i-im brightness profiles (namely the brightness versus the dis- 
tance from the center of the envelope) and SED with the observa- 
tions (described in the previous section) allows to constrain the 
density profile and, consequently, the temperature profile of the 
envelope. To be compared with the observations, the theoretical 
emission is convolved with the beam pattern of the telescope. 
Following the recommendations for the relevant telescope, the 
beam is assumed to be a combination of Gaussian curves: at 850 
yum, we use HPBWs of 14.5", 60", and 120", with amplitudes 
of 0.976, 0.022, and 0.002 respectively; at 450 yum, the HPBWs 
are 8", 30", and 12 0" with amplitude ratios of 0.934, 0.06, and 
0.006, respectively (ISandeU & Weintraubll200T!) : at 350 /urn, we 
use HPBWs of 9" and 22", with amphtude ratios of 0.8, 0.2, 
respectively (Section 2.1). 

In this work, we consider two cases for the density distri- 
bution. In the first one, we assum a broken power-law density 
profile as in the [Shul(ll977h solution: 

n(r) = n(ro) 1 — 1 r < n„f (1) 

n(r) = n(ro) j r > r,„/ (2) 

where r,„/ is the radius of the collapsing envelope (at larger 
radii the envelope is static). In the density profile, it repre- 
sents the radius at which the change of i ndex occurs and it is 
a free parameter Note that the complete IShul (Il977h solution 
also contains a transition part at the interface of the collapsing 
and static regions, just inside r,„/. In this region the slope gradu- 
ally changes from r ' to the limiting value of r"'^. A posteriori, 
the effective part of the envelope in r ' is relatively small and 
located in the inner region (< 1500 AU, equivalent to ~ 10")- 
Since the continuum observations, with resolutions of ~ 8-18", 
are not sensitive enough, we adopted the simplified structure de- 
scribed by the Eq. lll2l In the second case, we considered a single 
power-law density profile, where the index a is a free parameter: 

n(r) ^ n(ro){^J (3) 

In both cases, n(ro) is the density at ro, and the envelope starts at 
a radius /-,„ and extends up to ro„t- In total, both models have four 
free parameters determined by the best fit with the observational 
data: r,„y or a, n(ro), rou and r,„. Finally, DUSTY requires the 
temperature of the central source, T,, here assumed to be 5000 
K. Note that we verified that the choice of this parameter does 
not influence the resu lts, as already noticed by other authors (e.g. 
|j0rgensen et al.ll2002l) . In practice, the DUSTY input parameters 



Parameter 


Range 


Step 


a 


0.2-2.5 


0.1 


Y" 


50-2000 


10 


I'ln/ 


5-200 


2 




0.1-10. 


0.1 


T,„ 


300 K 


Fixed 


T. 


5000 K 


Fixed 



Table 2. Range and step of the DUSTY input parameters. 



are the infall radius r,„/ or the power-law index a, the optical 
thickness at 100 /zm, tioo, the ratio between the inner and outer 
radius, Y (=ro„,/r/„) and the temperature at the inner radius TuQ. 
The optical thickness is, in turn, proportional to n(ro) and r„,„. In 
both models, we obtain a lower limit to T,„ of 300 K, any larger 
value giving similar results. 

In addition to the above parameters, the opacity of 
the dust as function of the wavelength is a hidden pa- 
rameter of DUSTY. Following numerou s previous s tudies 
( Van Der Tak et al. 1999; E vans et al.ll2()()lt IShirlev et al.il2002l: 
Young et al.. .2003.: .Schoier et al. 11200^ I), we adop ted the dust 
opacity calculated bylOssenkopf & H ennind (1 19941) . specifically 
their OH5 dust model, which refers to grains coated by ice. 
Again, the basic result, though, does not substantially depend 
on the choice of the dust opacity model. 

As explained in llvezic & EhtzuJ (Il997h . DUSTY gives 
scale-free results, so that the source bolometric luminosity Lio/ 
and the distance are required to compare the DUSTY output 
with actual observations. We assumed here the lastest estimate 
of the distance to the p Ophiuchus cloud, namely 120 pc (see 
Introduction), and we derived the bolometric luminosity by inte- 
grating the observed emission over the full spectrum, optimizing 
the resulting (see below). 

We run grids of models for both cases described above. The 
summary of the covered parameter space is reported in Table |2] 
The best fit model has been found minimizing the with an it- 
erated two-steps procedure (see also Crimier et al. 2009). First, 
we use the observed brightness profiles at 350, 450 and 850 yum 
to constrain Y and a (or r,„y in the case of the Shu-like model), 
assuming a value for tioo- Second, we constrain the optical thick- 
ness Tioo by comparing the computed and observed SED, assum- 
ing the a (or r,„/) and Y of the previous step. The new tioo is 
used for a new iteration and so on. In practice, the iteration con- 
verges in two steps. This is because the normalized brightness 
profiles depend very weakly on tioo, while they very much de- 
pend on the assumed size of the envelope and on the slope of the 
density profile (see also the discussion in |j0rgensen et al.ll2002l: 
ICrimier et al.ll2009l) . In contrast, the dust optical thickness de- 
pends mostly on the absolute density of the envelope. Note that 
we constrain also the bolometric luminosity based on the best fit 
of the SEDB 

2.3. Results 

a) Brightness profiles and SED analysis 



The temperature at the inner radius T,„ in fact defines the radius at 
which the integration starts. The DUSTY codes requires the tempera- 
ture rather than the radius because it is based on a scale-free algorithm. 

' Note that the dust optical thickness Tioo affects the shape of the 
SED, so in general it enters in the determination of the bolometric lu- 
minosity. 
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1.3mm 


3mm 


Source 


F,(Jy) FWHM(") 


F,(Jy) FWHM(") 


A 
B 


0.49 1.18 
0.97 0.86 


0.16 1.62 
0.25 0.78 



Table 4. Integrated flux and Full Width at Half Maximum 
(FWHM) of the gaussians used to model the emission of the 
sources A and B at 1.3 and 3 mm. 



Table |3] lists the best fit parameters for the two models de- 
sc ribed in the previous s ection as well as the parameters obtained 
bv lSchoier et al.l (l2002l) for comparison. Figure [T] shows the so- 
lutions compared to the observed data (maps and SED). 

The results of this part of the analysis are: 

- The bolometric luminosity, obtained minimizing the is 
22 Lq rather then 14 Lq, the value obtained by simply scal- 
ing the bolometric luminosity adopted by ICeccarelli et alj 
(f2000b) and Schoier et al. (2002) to a distance of 120 pc in- 
stead of 160 pc. 

- The single power-law and Shu-like density distributions can 
both reproduce the observations, the maps and the SED data, 
included the IR part of the spectrum (Figure [T]i. They give 
similar best fit values, although the Shu-like distribution 
fits better the 450 //m profile. Also the derived physical pa- 
rameters are substantially similar for the two models. 

- The observed data, including the Spitzer data between 20 
and 40 iim, do not necessitate the presence of a large cavity 
(iJ0rgensen et al. 2005). Our solution with a larger luminosity 
(22 instead of 14 L©) can reproduce as well the MIR Spitzer 
observations. Figure |2] shows the MIR Spitzer and IRAS ob- 
servatio ns against the emission pre dicted by our models and 
these bv lj0rgensen et all (l2005h and lSchoier et alj (|2002|) . 

Figure [3] shows the derived dust temperature and density pro- 
files. The single power-law index density distribution predicts a 
slightly warmer and denser region at radii lower than about 100 
AU, but the differences are relatively small (see also Table[3]l. 

h) Interferometric data analysis 

The previous analysis considers single dish data, which at 
best have a spatial resolution of 8", equivalent to a radius of 
~ 1000 AU. In order to constrain better the inner region, we 
used 1 and 3 mm continuum interferometric observations ob- 
tained with the Plateau de Bure interferometer (PdBI) described 
in the previous section. In order to compare the model predic- 
tions with the observations, we have produced synthetic maps 
for each model, in which we added the predicted emission from 
the envelope plus the emission of the two Gaussian-like sources 
A and B. The parameters used to represent the emission of the 
sources A and B are extracted from the interferometric maps 
at 1.3 and 3 mm and reported in Table H] Then we used the 
UV_FMODEL task in Gilda^ to compute visibility tables with 
the same uv plane coverage than the actual observations. 

Table |5] reports the x^^^j values obtained using the observed 
and modeled visibility amplitudes at 1.3 mm. The comparison is 
done for each of the three models of the envelope, i.e. the single 
power-law density distributi on, the Shu - Uke de nsity distribution, 
and the central cavity of , Schoier et al.l (l2004l) . We tested cases 
in which the envelope is centered on one of the two sources or 
on the mid-way point between the two sources. The x^^d values 
are computed over more than lO** points. The 2D representation 

^ http://iram.fr/lRAMFR/GILDAS/doc/html/mis-html/nodel3.html 



of the visibility amplitudes observed and modeled in uv plan is 
difficult to read. Therefore to illustrate the comparison between 
the visibility amplitudes observed and modeled, we averaged the 
visibility amplitudes over the same uv radius and plotted them in 
Figure |4] Note that the error bars shown in Fig.|4]represent only 
the measurement eiTors. In fact, the standard deviations result- 
ing from the radial average are meaningless in this case because 
of the non-sphericial geometry (caused by the presence of the 
sources A and B). Figure |4] only aims to illustrate the compari- 
son between the observations and models, our conclusions will 
be based on the x^.^^ in Table |5] The computations were done 
considering all the visibility amplitudes obtained in the uv plane 
before the azimuthal average. 

The interferometric data are dominated by the two compo- 
nents of the IRAS 16293 binary system. However, regarding the 
envelope contribution — which is relevant for the present work — 
in the case where the envelope is centered on the mid-way point 
between the sources A and B, the visibilities at uv radii lower 
than about 80 kA are not well reproduced by the single power- 
law density profiles. Our solution with the Shu-like density 
profile and the solution with a cavity (suggested by Schoier et al] 
2004) give a much better agreement. However, when envelope 
is centered on B or A all three models give similar results (Fig. 
HJ, and fit well the observed visibilities. Since with the data 
available in the literature there is no way to know whether the 
envelope is centered on one of the two sources (A or B) or just 
on the mid-way point between them, either of the two solutions 
(centered on a source or on the mid-way point between A and 
B) is equally plausible. In other words, we are in a degenerated 
case dominated by the parameters relative to the binarity of the 
source. The interferometric data can be reproduced without a 
cavity. Note that the 3 mm PdBI observations do not provide 
additional information. We emphasize that the above analysis of 
the interf erometric data is d irectly applicable to the OVRO data 
that led Schoier et al. (2004) to suggest the presence of a cavity. 
The OVRO data probe visibility amplitudes lower than about 60 
k/l, a range also probed by the PDBI data. Therefore, as shown 
above, they can be reproduced by assuming that the envelope 
is centered on either A or B source without the necessity of 
assuming a cavity. 

c) Summary of the continuum data analysis 

Both the single dish and interferometric continuum data can 
be reproduced by either of the two models we considered: a 
Shu-like and single power-law density distribution. Furthermore, 
no cavity i s required to explain the data, neither the Spitzer 
MIR data (iJcirgensen et al.l l2005h nor the interferometric data 
(ISchoier et al.ll2004l) . 

Therefore, in the following we will adopt the Shu-like den- 
sity profile (which has a physical interpretation), with no cavity, 
as our reference model for the study of the gas temperature and 
water line predictions. 

3. Gas temperature profile 

3. 1 . Adopted method 

We compute d the gas tempera t ure pr ofile using the CHT96 code 
described in 'Ceccarell i et alJ (Il996l) (see als o ICeccarelli et al.l 
2000a; Maret et al. 2002t ICrimier et al.ll2009l) . Briefly, die code 
computes the gas equilibrium temperature at each point of 
the envelope, by equating the heating and cooling terms at 
ea ch point of the envelo pe. Following the method described 
in ICeccarelli et al.l (Il996i) . we considered heating from the gas 
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Quantity 


Shu-like model 


Power-law model 


ocnoier et ai. [ zuuz ) 




Best fit parameters 






22 


22 


on 

1.1 


D(pc) 


120 


120 


160 


TlOO 


2.0 


3.0 


4.5 


r,„f (AU) 


1280 






a 




1.8 


1.7 


Y 


280 


260 


250 


Physical quantities 


r.n (AU) 


22 


27 


32 


ro,, (AU) 


6100 


6900 


8000 


r(T<i„„=100K) (AU) 


76 


85 




n(Trf,,„=100K)(cm-3) 


2 X 10** 


3 X 10** 




M„„, (Mo) 


1.9 


2.1 


5.4 



Table 3. Summary of the dust radiative transfer analysis. The upper half of the t able lists the best fit parameters, the lower half 
reports some relevant physical quantities. The last column reports the results of the lSchoier et al.l (|2002|) analysis for comparison. 



10,00 



1,00 

0.10 

0,01 

15 20 25 30 35 40 
A (/i.m) 

Fig. 2. MIR Spitzer and IRAS fluxes as a function of the wavelength. The two crosses represent the integrated fluxes at 23.7 and 
25 jum obtained with MIPS aboard Spitzer and CPC aboard IRAS, respectively. The solid thick line is the IRS SPITZER spectrum 
observed between ~ 21-37 /vm (see text for details). The solid thin line is the best fit of the maps and SED data given by the model 
of envelope with a density profile with a broken power law with a'o«/=2 and a,,, = 1.5 in the outer and inner envelope respectively 
(Shu-like density distribution). The dashed line is the best fit given by the model of envelope with a single po wer-law index, a- \.%. 
The dotted-dashed l ine and the dotted lines show the emission resulting from the models obtained bv ,J0rgensen et al.l ( l2005h and 
ISchoier et alj (|2002|) . respectively. 




compression (due to the collapse), dust-gas collisions and photo- 
pumping of H2O and CO molecules by the IR photons emitted 
by the warm dust close to the centefl The cooling terms are 
the line emission from H2O, CO and O. The dust-gas collisions 
are a source of heating when the dust temperature is higher than 
the gas temperature and a source of cooling elsewhere. To com- 
pute the cooling from the li nes we used the code d escribed in 
ICeccarelU et ail(ll996ll2003h and lParise et aP (l2005h . The same 
code has been used in several past studies, whose results have 
been substantially confirmed by other groups (e.g. the analysis 
on IRAS 16293-2422 bv lSchoier et al- 2002). Briefly, the code is 
based on the escape probability formalism in presence of warm 
dust (see iTakahashi et al]|1983h . where the escape probability y6 



' Cosmic rays ionization is a minor heating term in the protostellar 
envelopes. 



is computed at each point by integrating the line and dust ab- 
sorption over the solid angle O as follows: 



fed 



+ fed (^L + fcd)* 



- exp [- (^L + ^d) ALth] 
ALth 



(4) 



where k-^ and fed are the line and dust absorption coefficients re- 
spectively, and ALth is the line trapping region, given by the fol- 
lowing expressions: 



ALth = 2Avi 



(5) 



in the infalling region of the envelope (where arcos \i is the angle 
with the radial outward direction) and in the static region (where 
7?env is the envelope radius): 



ALtf 



(6) 
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Fig. 3. Dust temperature (upper panel) and H2 density (lower panel) profiles from the best fit obtained in the two cases Shu-like 
(solid line) and single power-law density distribution (dashed line). 



In addition, H2O and CO molecules can be pumped by absorp- 
tion of the NIR photons emitted by the innermost warm dust. 
Since the densities and temperatures of the regions of the enve- 
lope targeted by this study are not enough to populate the levels 
in the vibrational states, the effect of the NIR photons is an extra 
heating of the gas, as described in Ceccarelli et al. (1996). Note 
that the code takes into account the dust with temperatures up to 
1500 K. 

The code has a number of parameters, that influence the gas 
equilibrium temperature. First, the dust temperature, assumed to 
be the output of the previous step analysis (®. Second, since 
the infall heating depends on the velocity gradient across the en- 
velope and the cooling depends on the gas lines (which can and 



are sometimes optically thick), the velocity field across the en- 
velope is also a parameter of the code. Here we assumed that 
the infall velocity field corresponds to the free fall velocity field 
with a 2 M0 source at the center of the envelope (ILoinard et al.l 
I2OO9I) . Third we adopted the "standard" cosmic rays ionization 
rate, namely 3 x 10"'^ s"', however, in practcie this parameter 
is unimportant. Finally, given the contributions of the H2O, CO 
and 01 to the gas cooling, their respective abundances are impor- 
tant parameters of the model. Previous theoretical studies have 
shown that the O abundance is constant across the envelope, ex- 
cept in the very inner regions, where the gas temperature exceeds 
about 250 K and endot hermic reactions that form O H and H2O 
become very efficient (iCeccarelli et al.1 [1996; ,Dotv & Neufeldl 
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Fig. 4. Visibility amplitudes as a function of the uv radius of the 1.3 mm PdBI observations (squares) against the predictions 
ob tained with the best fit models of a Shu-like (solid lines) and single power-law (dashed lines) density distribution. The model 
of ISchoier et alJ (l2004l) . with the central cavity, is shown for a comparison (dotted lines). Panels (a) and (b) show the case where 
the envelope is centered on the mid-way point between the sources A and B. Panels (c) and (d) show the case where the envelope 
is centered on the source B. Panels (a) and (c) show the visibility amplitudes modeled with the synthetic maps containing the 
contribution of the envelope and the sources A and B. Panels (b) and (d) show the visibility amplitudes modeled with the synthetic 
maps containing the contribution of the envelope only. Finally, the zero-signal expectation values are represented by the empty 
triangles on the c panel. 



[T997HDotv et al.l I2OO4I) . Unfortunately, this parameter is very 
poorly constrained by observations (because of the difficulty of 
observing the O fine structure lines and the fact that they are 
easily excited in the foreground molecular cloud and its asso- 
ciated PDR (see e.g. Liseau et al.lll999HCaux et al Here 
we assumed that the atomic O abundance is equal to 1 x 10"^ 
and verified a posteriori that the O fine structure line emission is 
consistent with the ISO observations. We note, however, that O 
is never the dominant coola nt except, perhaps, in a very small re - 
gion of the envelope ( Ceccarelli et al.i 2000aV Maret et ai]|2002h . 
CO, in contrast, is the main coolant in the outer envelope, but, 
since the cooling lines are heavily optically thick, its abundance 
does not play an important role (in the regime where it is higher 



than about 1x10^ with respect to H2). We, therefore, assumed 
a canonical 1 x 10 CO abundance across the envelopeQ 

Finally, water is an important coolant in the inner region, 
where the grain mantles sublimate, injecting into the gas phase 
large quantities of water molecules. Again, based on previous 
studies, we approximated the H2O abundance with a step func- 
tion: the abundance is X,„ in the region where the dust temper- 
ature exceeds 100 K, and Xo„, elsewhere. Both X,„ and Xo„, are 
found by comparing the theoretical predictions with the ISO ob- 
servations (see ^3.2l i. Note that, to solve the water level popula- 
tion statistical equilibrium equations, we used the collisional co- 



Note that it possible that a region where CO abundance is lower, 
because of the freezing onto the dust grains, exists. However, the gas 
cooling is relatively insensitive to this lower abundance for the reasons 
explained in the text. 
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Fig. 5. Reduced-;\f^ contours obtained comparing the model pre- 
dictions and the observations towards IRAS 16293, as function 
of the inner X,„ and outer X,,,,, H2O abundance. In this compu- 
tations, we used the dust temperature and density profiles of the 
Shu-like model of Table [3] 
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Fig. 6. Ratio of the observed over predicted H2O line flux versus 
the upper level energy in cm"'. In this computations, we used 
the dust temperature and density profiles of the Shu-like model 
of Table [3] and an H2O abundance equal to 5 x 10"^ in the inner 
envelope and 2.5 x 10"^ in the outer envelope (see text). 



efficie nts between H2O and H2 recently computed bv lFaure et alJ 
(I2007h . We assumed that the ortho-to-para H2 ratio is at the 
Local Thermal Equilibrium (LTE) in each part of the envelope. 
Finally, we assumed a H2O ortho-to-para ratio equal to 3. 

3.2. Results 

In order to constrain the water abundance, necessary to predict 
the gas temperature profile, we used the ISO observations of the 
water lines towards IRAS16293 (Ceccarefli et al. 1998, 2000a) 
and compared with the model predictions obtained with differ- 
ent H2O abundance profiles. For that, we run a grid of models 
with X,„ and Xout varying between 1 x 10 and 1 x 10"^. Figure 
|5] shows the resulting as function of the inner and outer H2O 
abundance, in the case of the Shu-like structure of Table [3] The 
H2O abundance in the outer envelope is constrained between 
0.7-20 xlO"^ at 3 cr, with the best fit solution at 2.5 x 10"^ 
(reduced 1.3). The inner H2O abundance is even less con- 
strained, lower than about 2 x 10"^ at 3 cr confidence level, and 
0.5x10"^ at 1.5 cr confidence level, i.e. slightly lower than the in- 
ner abundance, which is in contradiction with the hypothesis that 
ice mantles sublimate injecting water in the gas phase. However, 
one has to consider that the inner H2O abundance is very poorly 
constrained by the ISO observed lines, which are very optically 
thick and not high enough in energy (as it is shown by the pre- 
dictions of Table|7]i. We, therefore, only consider significant the 
3cr level limit to the water inner abundance, warning that even 
that may be questionable for the same reasons (see also ® . The 
ratio between the observed by ISO and predicted lines fluxes as 
function of the upper level energy is shown in Fig.|6] 

The results for the single power-law density profile are simi- 
lar, with the upper limit on the inner H2O abundance a factor ten 
lower (because of the higher density in the inner part predicted 
by this model). 

Figure |7] shows the dust and gas temperature profiles for the 
Shu-like density distribution. Gas and dust are thermally coupled 
across the whole envelope with the largest difference (< 10%) in 
the infall region, where r < r,„/. Figure [8] shows the heating and 
cooling terms across the envelope. The heating is totally domi- 
nated by the compression in the entire collapsing region, and by 



the collisions between the dust and the gas in the outer envelope. 
Note that the H2O and CO photo-pumping plays only a minor 
role. The cooling is dominated by the water line emission in the 
inner region (where the ice sublimate), by the dust-gas collisions 
in a large intermediate region, and by the CO line emission in the 
very outer envelope. These results are very similar to those of 
fCeccarelli et al. (2000a). Note that low-lying water lines could 
possibly be contaminated by the outflow driven by the source 
This effect could possibly lead to an overestimation of the outer 
water abundance. Since the ISO water lines are spectrally unre- 
solved, it is not possible to address this question using the ISO 
data. 

Tables |7] and |6] list the predicted fluxes of the water lines 
which will be observable with HSO, for the best fit model (H2O 
abundance in the outer envelope equal to 2.5 x 10"^) with H2O 
abundance in the inner envelope equal to 0.5 x 10"^, but also for 
the case with a larger H2O inner envelope abundance (2 x 10"^). 

4. Discussion and conclusions 

The new analysis of the single dish and interferometric contin- 
uum observations of the envelope of IRAS 16293 confirms that 
an envelope of about 2 Mq surrounds the proto-binary system 
of IRAS 16293. The envelope can be described with a Shu-like 
density distribution, corresponding to the gas collapsing towards 
a 2 M0 central star. The luminosity of IRAS 16293 has been re- 
evaluated to be 22 L0 for a distance of 120 pc. 

Both the single dish and interferometric data can be repro- 
duced by an envelope with an inner radius between 20 and 30 
AU, equivalent to about 0.4", and smaller than the radius at 
which the dust temperature reaches 100 K (the ice mantle sub- 
limation temperature, namely 75 to 85 AU). We found that our 
new analysis can reproduce the full SED, including the Spitzer 
MIR data, without the nece ssity of a central cavity of 800 AU ra- 
dius (|j0rgensen et alJl2005h . The diff'e rence between our models 
and the previous ones (based on the ISchoier et al.i r2002 initial 
model) is the larger bolometric luminosity (kept as a free pa- 
rameter in our models) and the lower optical thickness of the 
envelope ( tioq^,,,-!, namely twice smaller than in the Schoier's 
models). These differences contribute to make that the predicted 
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Fig. 7. Dust (solid line) and gas (dashed line) temperature profiles of the best-fit obtained assuming a Shu-like density distribution 
(see text). 



MIR radiation flux agrees with the observed one. Finally, the in- 
terferometric data, being dominated by the two components of 
the binary system, do not provide significant constraints on the 
envelope structure, with one exception. They exclude the case of 
an envelope with a single power-law density profile centered on 
the mid- way point between th e two sources. Based on this as- 
sumption, ISchoier et al.) (l2004 suggested the presence of a cav- 
ity 800 AU in diameter. Since no data constrain where the center 
of the envelope is located, we favor the solution with the enve- 
lope centered on one of the two sources. In addition, the Shu-like 
model fits slightly better the PDBI data and is thus our preferred 
solution. 

As already noted by ICeccarelli et al.1 (l2000al) . the ISO data 
do not allow to constrain the inner H2O abundance because the 
detected lines are optically thick and cover a relatively low range 
of upper-level energies (< 300 cm"'). Also the outer-envelope 
water abundance is relatively poorly constrained. In addition, 
the relatively low spectral resolution of ISO does not allow to 
determine whether some lines are contaminated by the emission 
from the outflow. The future observations with the HIFI spec- 
trometer aboard the Herschel Space Observatory, launched in 
May 2009, will certainly constrain better the water abundance 
profile across the IRAS 16293 envelope, helping to understand 
the distribution of water in protostars similar to the Sun progen- 
itor. On the one hand, the water abundance in the outer enve- 
lope (0.7-20x10"^) derived from the ISO observations here is 
consistent with some previous estimates of water abundance in 
cold gas (e.g. Cernicharo et al. 1997) but only marginally consis- 
tent with other low estimates (e.g. lSnell et al.l l2000). so the new 
Herschel/HIFI observations, with their largely improved spatial 
and spectral resolution, will be crucial in settling the question. 
On the other hand, the inner envelope abundance (<2xl0"^) 
is lower than expected if all ice in the mantles sublimates (see 
ICeccarelli eta l. 2000a). Also in this case, the Herschel/HIFI ob- 
servations will help to understand this point. 



Finally, as stated in the Introduction, the major scope of the 
present work is to provide as accurate as possible estimates of 
the dust and gas temperature profiles of the cold envelope of 
IRAS 16293 and its warm inner component, also known as the 
hot corino, to interpret the data observed in two large projects, 
TIMASSS and CHESS (see Introduction). We are aware that the 
proposed description has the intrinsic and clear limit of not tak- 
ing into account the multiple nature of the IRAS 16293 system. 
But it has the merit in allowing the interpretation of the single- 
dish observations of the upcoming projects, within this limita- 
tion. A more detailed analysis will only be possible once the rel- 
evant molecular emission is observed with interferometers, re- 
solving the two components of the system. In absence of that, 
the analysis based on a single warm component and cold enve- 
lope is the only vialable and allows a first understanding of the 
chemical composition of a system which eventually will form a 
star and planetary system like our own. 
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in between A and B 


onB 


on A 


Model 


2 




'I 

X ml 


Single power-law 


18.5 


10.2 


8.5 


Shu-like 


13.3 


10.3 


8.3 


Central cavity 


11.0 


11.3 


9.9 



Table 5. x^^^ values obtained comparing the visibility amplitudes observed and modeled at 1.3 mm. The number of degrees of 
freedom v is 10230. The 1" column reports the;^^^^^^ obtained using the synthetic maps with the envelope centered on the mid-way 
point between the sources A and B. The 2"'' and the y' columns report ihsxled obtained using the synthetic maps with the envelope 
centered on the source B and on the source A, respectively. The 1 the 2""', and the y' lines report thtx^ed obtained in all the cases 
discussed in the text, using the model o f envelope with a single power-law density distribution, a Shu-like density distribution, and 
a central cavity of lSchoier et al i (|2004|) . respectively. 
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HIFI range 




^in{^2^) — ^ X iU 


■'^in(,n2(J;— J X iU 


Transition 


Frequency 


Flux 


Flux 






(Jv Km S ; 


(K. Km s ) 


lio — ^ loi 


557.0 


32 


32 


532 441 


020.7 


0.31 





2ii 2q2 


752.0 


23 


23 


422 hi 


916.2 


0.04 


0.04 


2o2 — > 111 


987.9 


37 


37 


3l2 — > 3o3 


1097.3 


26 


26 


111 — > Ooo 


1113.4 


46 


46 


^25 ~> 818 


1 146.6 


0.04 





3l2 221 


1 153.1 


34 


34 


634 541 


1158.3 


1.0 





321 — » 3l2 


1162.9 


13 


11 


422 4i3 


1207.6 


1.1 


1.1 


220 — > 2n 


1228.8 


7.6 


7.6 


523 ~> 5i4 


1410.7 


5.7 


2.1 


643 — > 7i6 


1574.2 


0.10 





4l3 ^ 4o4 


1602.2 


2.2 


2.2 


221 — > 2i2 


1661.0 


40 


37 


2l2 ~> loi 


1669.9 


66 


67 


432 5o5 


1713.9 


2.1 


0.1 


3o3 ~* 2i2 


1716.8 


56 


56 


'34 ^ '25 


1 7Q7 




n 1 

U. i 


532 — * 523 


1867.7 


4.6 


0.9 


634 — > 7o7 


1880.8 


0.41 





845 ~* 752 


1884.9 


0.09 





331 4o4 


1893.7 


0.02 


0.02 



Table 6. Predictions of the line fluxes (after subtraction of the continuum) of the water lines observable with the Herschel spectrom- 
eter HIFI. The third (fourth) column reports the predictions computed assuming that the inner H2O abundance is equal to 2 X 10~* 
(5 X 10 **) and the outer H2O abundance is 2.5 x 10"^, for the Shu-like model best fit. 
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PACS range X,„(H20)^2 x 10" X,„(H20)^5 x 



Transition 


Wavelength 


Flux 


Flux 


Flux 






(lO-'^erg cm-2) 


(lO-'^erg cm-2) 


(lO-'^erg s-i ( 


936 — > 845 


62.42 


0.05 







9l8 9o9 


62.93 


0.4 







818 ^ 7o7 


63.32 


1.0 


0.1 




661 652 


63.91 


0.2 







761 — > 752 


63.96 


0.05 







625 — > 5i4 


65.17 


1.7 


0.5 




7l6 ^ 625 


66.09 


1.2 


0.2 




330 221 


66.44 


2.5 


1.6 




331 220 


67.09 


0.3 


0.3 




330 — > 3o3 


67.27 


2.0 


0.5 




827 — > 818 


70.70 


0.6 







524 ^ 4i3 


71.07 


0.03 


0.03 




7o7 616 


71.95 


1.1 


0.2 




725 — > 634 


74.95 


0.6 







321 — > 2i2 


75.38 


2.9 


2.4 


2.7±1.0 


652 643 


75.83 


0.3 







550 — > 541 


75.91 


0.5 







752 743 


77.76 


0.2 







423 ^ 3i2 


78.74 


2.0 


1.5 




927 ^9,8 


81.41 


0.10 







616 5o5 


82.03 


1.1 


0.5 




836 — * 827 


82.98 


0.2 







606 ^ 5i5 


83.28 


0.04 


0.04 




7l6 — » 7o7 


84.77 


0.8 


0.1 




845 ^ 836 


85.77 


0.06 







322 — > 2ii 


89.99 


0.5 


0.5 


0.5±0.5 


643 ^ 634 


92.81 


0.3 







625 616 


94.64 


0.7 


0.2 




441 ^ 432 


94.71 


0.6 


0.1 




5l5 — » 4o4 


95.63 


0.3 


0.3 




541 — » 532 


98.49 


0.5 







5o5 ^ 4i4 


99.49 


1.3 


1.1 




5 14 — > 423 


100.91 


1.1 


0.8 


1.3±0.6 


220 — > 111 


100.98 


0.8 


0.8 




634 625 


104.09 


0.3 







221 ^ lio 


108.07 


2.2 


2.1 


1.7±0.6 


743 734 


112.51 


0.07 







4i4 — > 3o3 


113.54 


1.5 


1.3 




734 643 


116.78 


0.10 







432 423 


121.72 


0.4 


0.2 




4o4 — > 3i3 


125.36 


0.6 


0.6 




331 — > 322 


126.71 


0.1 


0.1 




725 ^ 7i6 


127.88 


0.10 







423 ^ 4i4 


132.41 


0.7 


0.6 


0.7±0.7 


836 — > 743 


133.55 


0.03 







5i4 ^ 5o5 


134.94 


0.5 


0.3 




330 — > 321 


136.49 


0.6 


0.4 




3i3 — » 2o2 


138.53 


0.9 


0.9 


0.9±0.7 


4l3 — » 322 


144.52 


0.1 


0.1 




322 — > 3i3 


156.20 


0.2 


0.2 




523 ^ 432 


156.26 


0.3 


0.04 




532 — > 523 


160.51 


0.1 







734 I25 


166.81 


0.05 







3o3 ^ 2i2 


174.62 


1.2 


1.2 


2.5±0.6 


432 5o5 


174.92 


0.05 







2l2 — > loi 


179.53 


1.4 


1.4 


2.9+0.6 


221 2i2 


180.49 


0.9 


0.8 


0.9+0.4 


4l3 ^ 4o4 


187.11 


0.05 


0.05 
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Table 7. Predictions of the line fluxes (after subtraction of the continuum) of the water lines observable with the Herschel spectrometer PACS. The 
third (fourth) column reports the predictions computed assuming that the inner H2O abundance is equal to 2 X 10~* (5 X 10"*) and the outer H2O 
abundance is 2.5 x 10 for the Shu-like model best fit. The fifth column reports the observed ISO fluxes 



